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Abstract 



The name Ising has come to stand not only for a specific model, but 
for an entire universality class - arguably the most important such class 
- in the theory of critical phenomena. I review several examples, both 
in and out of equilibrium, in which Ising universality appears or is 
pertinent. The "Ornstein-Zernike" connection concerns a thermody- 
namically self-consistent closure of the eponymous relation, which lies 
at the basis of the modern theory of liquids, as applied to the Ising 
lattice gas. Debye and Hiickel founded the statistical mechanics of 
ionic solutions, which, despite the long-range nature of the interaction, 
now appear to exhibit Ising-like criticality. The model of Widom and 
Rowlinson involves only excluded-volume interactions between unlike 
species, but again belongs to the Ising universality class. Far-from- 
equilibrium models of voting behavior, catalysis, and hysteresis provide 
further examples of this ubiquitous universality class. 
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I. INTRODUCTION 



Since the playful title of this review might generate confusion, let me start by 
saying that I have no idea whether or not Ernst Ising actually met any of the other 
scientists mentioned. The world lines intersecting here belong to models and theo- 
ries, not to persons, real or imagined! Much as the name Galileo has come to define 
a specific kind of relativity, Hamilton a kind of dynamics, and Gauss a distribu- 
tion, so "Ising" is now indelibly associated with a specific kind of critical behavior, 
a "universality class" in the familiar jargon of renormalization group theory. The 
latter tells us (as is borne out by experiment, analysis, and simulation), that the 
factors determining scaling properties in the neighborhood of the critical point are 
few, and pertain to very general properties of the system, such as its dimensionality 
and that of the order parameter, range and symmetry of interactions. For historical 
reasons, the class defined by a scalar order parameter, 0, short-range interactions, 
isotropy, and symmetry under inversion (<f) — > —<f) in the absence of an external 
field), is commonly known as the Ising universality class. While we should expect 
many systems to fall in this class, most of the present article is devoted to examples 
whose membership is surprising or controversial. 

My first example, self-consistent Ornstein-Zernike theory (Sec. II), is not a bona- 
fide member of the Ising class, but is, in a sense, making an impressive effort to be 
one! Criticality in electrolyte solutions, discussed in Sec. Ill, remains controversial, 
though recent theoretical and experimental studies support Ising-like behavior. We 
will see that the simplest lattice version of the well-known primitive model may 
exhibit a line of Ising-like critical points. We are used to thinking of the Ising model 
as possessing two relevant variables (or scaling fields), the temperature and the 
external magnetic field. The Widom-Rowlinson model (Sec. IV) (in its lattice-gas 
version) illustrates the possibility of Ising criticality in a system with no temperature 
to speak of. Having dispensed with temperature, we proceed in Sec. V to the more 
radical step of eliminating the Hamiltonian (and thermodynamics) altogether, in 
several far-from-equilibrium stochastic processes. The guiding principle of symmetry 
again permits one to understand the appearence of Ising criticality in these systems. 
A brief summary is provided in Sec. VI. 

II. SELF-CONSISTENT ORNSTEIN-ZERNIKE THEORY FOR THE 

ISING MODEL 

Ornstein-Zernike theory |1| represents an important, early chapter in the theory 
of critical phenomena, and is still widely used to interpret data from scattering 
experiments. The Ornstein-Zernike relation (OZR), moreover, serves as the starting 
point for modern theories of liquids . This relation introduces the direct correlation 
function c(r) |4j in terms of the total correlation function h(r) = g(r) — 1, where 
g(r) is the radial distribution function: 




(1) 
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Here p is the number-density. Of course one must now supply a closure, i.e., a second 
relation between c(r), the intermolecular potential w(r), and (in certain cases) h(r) 
itself. This, in general, can be done only in a heuristic manner, as for example in 
the mean-spherical approximation, c(r) = —w(r)/kBT, or the more sophisticated 
Percus-Yevick or hypernetted chain closures [0|§. Surprisingly, these simple closures 
lead in many cases to good predictions for h(r), and for thermodynamic properties, 
via one of the standard, exact relations (the so-called virial, compressibility, and 
energy routes) involving h(r). This is remarkable, since neither the OZR nor the 
closure include any thermodynamic input. The OZ formalism and the resulting 
'integral equations' for h(r) are couched purely in terms of correlations, without 
reference to a free energy. 

The absence of thermodynamic input becomes evident when one compares the 
predictions from different routes. Thus the Percus-Yevick closure yields two dif- 
ferent equations of state for the hard-sphere fluid, depending on whether the link 
to thermodynamics is made using the virial or the compressibility route. Such 
inconsistencies suggest a new strategy for closing the OZR, in which we use a ther- 
modynamic relation, expressing equality of (for example) the pressure evaluated via 
different routes, to find a second relation between h, c, and w. Such an approach 
was proposed over twenty years ago by H0ye and Stell |J, but detailed numerical 
implimentations have appeared only in the last few years. The first such application 
was, naturally, to the three-dimensional Ising model, which we outline here. (For 
details the interested reader may consult Refs. and 0.) 

We consider the Ising lattice gas with independent variables p (number of par- 
ticles per site) and (3 = 1/ksT. The interparticle potential w(r) is zero for r > 1, 
-1 for r = 1 (an attractive nearest-neighbor interaction) and infinite for r = 0. The 
latter implies the 'core-condition' h(0) = —1, expressing the simple fact that two 
particles may not occupy the same site. On a lattice the virial route, which involves 
the spatial derivative of the potential, does not exist. Consistency between the two 
remaining routes, energy and compressibility, is embodied in the relation 

d 2 ( P u) = d 2 ((3p) 

p dp 2 dpdp ' 1 ' 

where u is the internal energy per particle and p is the pressure. The energy relation 
is 

u = ~qp{l + hi), (3) 

where hi denotes the total correlation function at the nearest-neighbor separation, 
and q is the coordination number. The inverse compressibility is related to c via 

^ = l-*0) (4) 

where c denotes the Fourier transform. The direct correlation function is, in effect, 
defined by the OZR, which reads, on lattice, 
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hi = Ci + p^Ci-fa , 



(5) 



(the sum is over the sites of lattice C). This relation, or, more specifically, its Fourier 
representation, 

enables us to find h given c. Since we already have h(0) = — 1, the above relations, 
incorporating thermodynamic consistency, suffice for determining one free parame- 
ter, as a function of p and f3. At this stage we introduce the unique approximation 
of our theory, which is to set c(r) = for r > 1, i.e., beyond the range of the interac- 
tion. This defines what we call the self- consistent Ornstein-Zernike approximation 
or SCOZA. After some manipulations we obtain a nonlinear partial differential equa- 
tion for c(l) = ci(P,p) (c(0) is effectively fixed by the core condition). Integrating 
this PDE numerically eventually yields h(r) and all thermodynamic properties of 
the lattice gas. 

The results compare remarkably well against the best numerical (series expan- 
sion) estimates for the lattice gas [§,0. We find the critical temperature from the 
condition of a diverging compressibility at p = 1/2 (below T c , infinite compressibil- 
ity signals the spinodal); the coexistence curve is constructed in the usual manner. 
Especially impressive are the predictions for critical parameters, listed in Table I; 
SCOZA reproduces the best series estimates to within 0.2%. 
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TABLES 



TABLE I. Comparison of critical parameters 



lattice 


Pc,scoz 


fic,be 


UcySCOZ 


Uc,be 


Sc,scoz 


$c,be 


(Pp)c 


(Pp)c,be 


SC 


0.88497 


0.88662 a 


-2.010806 


-1.9961 6 


1.1037 


1.1158 6 


0.1140 


0.1124 


BCC 


0.62848 


0.62947 c 


-2.564460 


-2.5464 6 


1.1542 


1.1641 6 


0.1256 


0.1244 


FCC 


0.40772 


0.40825 c 


-3.768955 


-3.7423 6 


1.1703 


1.1804 6 


0.1302 


0.1292 



a Ref [8] b Ref [9] c Ref [10] 
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Can SCOZA predict critical exponents? Its thermodynamic predictions are so 
accurate that effective exponents (i.e., derivatives of the compressibility, order pa- 
rameter, etc., with respect to temperature) are in reasonable agreement with Ising 
model values near, but not asymptotically close to, the critical point. We really 
should not expect SCOZA to reproduce Ising critical exponents: c(r) is not zero for 
r > 1 in the Ising model. In particular, it develops a power-law tail at the critical 
point. The actual critical behavior of the SCOZA-lattice gas is unusual f|?J: for 
T > T c the exponents are those of the spherical model: 7 = 2, 5 = 5, and a = — 1 
(the specific heat does not diverge at T c ), as compared with 7 ~ 1.24, 5 ~ 4.8, 
and a ~ 0.1 for the three-dimensional Ising model. Below T c , SCOZA yields a new 
set of exponents: 7' = 7/5, (3' = 7/20 and a' = —1/10, considerably better than 
mean-field or spherical model values. 

Clearly the way to improve SCOZA is to relax the truncation of c(r). But even 
with its current limitations, SCOZA represents the first theory to yield globally 
accurate thermodynamic and structural properties for a fluid model. It has been 



applied with success to off-lattice fluids [11], and fluids in a disordered matrix [12 



and promisses to become an important tool in the study of liquids. 



III. ISING UNIVERSALITY AND THE PRIMITIVE MODEL OF 

ELECTROLYTES 

The critical behavior of electrolyte solutions remains a challenging area, both ex- 
perimentally and theoretically. Experiments on various ionic solutions have yielded 
either mean- field like critical exponents |13[ , Ising-like exponents [[14]], or a crossover 



from mean- field to Ising-like behavior as one nears the critical temperature | I5| . In 
some cases the crossover from mean-field to Ising exponents occurs at a reduced 
temperature much smaller than that typically seen in non-Coulomb liquids. 

Theoretical studies have tended to focus on the simple system introduced by 
Debye and Hiickel in 1923 [ 16(1 , commonly known as the restricted primitive model 



(RPM). Here the positive and negative ions are represented by hard spheres of diam- 
eter a, containing point charges ±g at their centers. The solvent is not considered 
explicitly; its effect is represented solely by the dielectric constant e that enters the 
expression for the electrostatic potential. While clearly a minimal model for ionic 
solutions, the RPM already presents great difficulty to theory and simulation. Debye 
and Huckel's analysis of the dilute, high-temperature regime led to the fundamen- 
tal result that electrostatic interactions are screened on a length scale P -1 , where 
T 2 = 4irpf3q 2 /e. (The effective interaction or potential of mean force takes on a 
Yukawa form oc e _r//r /r.) 

What sort of phase diagram should we expect for the RPM? The hard-sphere 
part of the potential alone will lead to a fluid-solid transition. At high temperatures 
(T* = aksTe/q 2 3> 1) we expect packing considerations to dominate, so that the 
high-density structure is FCC or HCP. But since these structures are incompatible 
with antiferromagnetic order (i.e., positive and negative charges occupying distinct 
sublattices), we should expect a structural phase transition to a bipartite lattice 
(presumably BCC) as we lower the temperature, at high density. 
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The opposite corner of the p — T plane, low temperature and density much 
less than the solid, is where we find the critical point in simple fluids. Since the 
RPM lacks the short-range attraction that drives the liquid-gas transition, it is not 
immediately obvious that it should exhibit such a transition. Nevertheless, Stell, 
Wu, and Larsen, using liquid-state theory, reached the conclusion that the RPM has 
a liquid-gas coexistence curve with a critical point [|17]], but with a critical density 
much lower than for a simple argon-like fluid. Monte Carlo simulations confirm 
these conclusions, but only in the last few years have the studies of different groups 
converged toward common values for the critical density and temperature of the 
RPM [p^-pO|. Given the difficulty in simply locating the critical point, it is not 
surprising that current simulations shed little light on its nature. 

On the theoretical side, however, there has been considerable discussion of RPM 



criticality, pointing toward Ising-like behavior |[2T|-[23|j. This seems at odds with 



conventional wisdom (Ising universality for short-range, and mean-field behavior for 
long-range interactions), but may be understood intuitively as follows. At the low 
temperatures of interest (note that T* « 0.05), the Coulombic interaction strongly 
suppresses charge-density fluctuations on scales larger than a; ions associate into 
pairs and larger aggregates. Thus the effective interacting units are not individual 
ions but clusters that are typically neutral (or nearly so), interacting via multipolar 
forces (presumably quadrupole and higher) that are of short range and, on average, 
attractive (since such fluctuations lower the energy). 

Such an intuitive picture finds support in recent experiments showing Ising crit- 
ical exponents at an ionic-solution critical point |T5|]. What is needed, from the 
theoretical standpoint, is an argument that takes us from the RPM to a continuum 
description of density fluctuations, Sp(r), having the same form (up to irrelevant 
terms) as that for the Ising model, i.e., the usual 4 field theory. It is clear, on 
the other hand, that the RPM needs to be described in terms of a pair of coupled 
fields, the mass or number density 5p and the charge density tp. An important first 
step is the proper formulation of a mean-field theory [24,25j. Very recently, Ciach 
and Stell constructed a Landau- Ginzburg free energy functional, starting from the 
mean-field theory of the RPM, in terms of the two fields, 5p and ip. They find 
that integrating out the charge density fluctuations leads to an effective field theory 
for Sp having the expected 4 form. (Essentially, an attractive effective interaction 
between mass-density fluctuations is mediated by charge-charge correlations.) Thus 
the crucial link between the RPM and Ising-like criticality appears to be at hand. 

In light of the great difficulty of RPM simulations, it seems useful to study 
a lattice restricted primitive model (LRPM) p^ , p7|l , since lattice simulations offer 
substantial economies in the computation time required for evaluating the poten- 
tial. (Overlap checks are trivial using a site-occupancy matrix, while the Coulomb 
potential - with the contributions from the infinite periodic array of cells suitably 
accounted for - may be stored in a lookup table.) Dickman and Stell considered 
a lattice gas of particles interacting via site exclusion (multiple occupancy forbid- 
den) and a Coulomb interaction u(rij) = SiSj/rij, where = |rj — r,-| is the distance 
separating the particles (located at lattice sites and r 3 -), and Si = +1 or — 1 is the 
charge of particle i. (Exactly half the particles are positively charged, half negative. 
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They are restricted to a simple cubic lattice with periodic boundaries. ) 

At full occupancy, the LRPM may be viewed as an antiferromagnetic Ising model 
with long-range interactions. One naturally expects a critical or Neel point sepa- 
rating a high-temperature phase from a phase exhibiting antiferromagnetic order. 
(The Ising analogy facilitates formulation of a mean-field theory for the LRPM.) 
Indeed, it was proven some time ago that the LRPM on the simple cubic lattice 
exhibits long-range order at sufficiently low temperatures and high fugacities [EB 



In simulations, we find that this point marks one terminus of a line of critical points 
in the p — T plane (see Fig. 1). The other end of the critical line is a tricritical 
point, which intersects the coexistence curve between a low-density disordered phase 
and the ordered phase. The Monte Carlo simulations used a relatively small lattice 
(16 3 sites), sufficient to map out the phase diagram by studying the order param- 
eter, specific heat, and correlation functions, but too small to yield information on 
critical behavior ||26|| . 

Thus the simplest version of the LRPM shows a phase diagram that is rather 
different than the continuous-space RPM, exhibiting a critical line and tricritical 
point rather than a simple critical point. One reason for this difference is that the 
simple cubic lattice facilitates antiferromagnetic order. It may be that on a lattice 
that frustrates such order (the FCC structure, for example), the critical line will 
disappear, and the tricritical point become a critical point. On the other hand, 
Panagiotopoulos and Kumar found exactly this, on the simple cubic lattice, when 
the hard-core exclusion range is > 3 times the lattice spacing |27|| . A suitably 
modified version of the LRPM may therefore be useful for studying ionic criticality. 
The phase diagram appears to be highly sensitive to changes in lattice structure or 
short-range interactions, as was also noted by Ciach and Stell . 

Even if the critical line of the simple-cubic LRPM is not exactly what we had 
in mind for understanding the off-lattice RPM, its nature, and that of the associ- 
ated tricritical point, are of interest, and potentially relevant to other systems with 
Coulombic interactions. Given the symmetry of the system, the Ising universality 
class again seems the natural candidate. The (very limited) simulation data seem 
to indicate (3 < 0.326, (the 3-d Ising value), which should in any case rule out a 
mean-field type transition. Clearly larger-scale simulations and finite-size scaling 
analysis are in order; a low-temperature expansion for the fully occupied case might 
also prove useful. We close this section with the observation that Ising criticality 
is compatible with a long-range bare interaction, provided the effective interactions 
between critical fluctuations are of short range. 



IV. ISING WITHOUT TEMPERATURE: THE WIDOM-ROWLINSON 

MODEL 

In this section we again consider a model imported to the lattice from its original 
continuous-space formulation. The Widom-Rowlinson (WR) hard-sphere mixture is 
perhaps the simplest binary fluid model showing a continuous unmixing transition, 
and has been the subject of considerable study regarding its thermal and interfacial 
properties P5| , |5DB , as has the Gaussian /-function version of the model introduced 



S 



somewhat earlier by Helfand and Stillinger [[H]]. Despite this interest, however, 
definitive results on the location and nature of the WR critical point are lacking. As 
a first step in this direction, Dickman and Stell performed extensive simulations of 
the lattice-gas analog of the WR hard-sphere mixture the Widom-Rowlinson lattice 
model (WRL) @. 

In the original WR model, AB pairs interact via a hard-sphere potential whilst 
AA and BB pairs are noninteracting |29[| . The WRL model is a two-component 
lattice gas in which sites may be at most singly occupied, and in which nearest- 
neighbor A-B pairs are forbidden. Like the WR model, this is evidently an athermal 
model (all allowed configurations are of the same energy), and is characterized solely 
by the densities of the two species or by the corresponding chemical potentials pa 
and /I b- {pa = Pb = Pc of course, at the critical point.) 

Despite the absence of a temperature or energy scale, the WRL model has a close 
affinity to the Ising model. To see this, note that the WRL may be viewed as an 
extreme member of a family of binary alloy models with nearest-neighbor interactions 
that are repulsive between unlike species (interaction energies €ab > 0, eaa = £bb = 
0). The Ising model may be transcribed into such a model by identifying up and 
down spins with A and B particles, respectively, yielding a "close-packed" alloy 
that unmixes at the Ising critical temperature. Allowing a small fraction of vacant 
sites results in a dilute binary alloy (DBA) with a somewhat depressed critical 
temperature; continuing the dilution process, one arrives at a model with T c = 0. 
This zero-temperature terminus of the DBA critical line is precisely the WRL critical 
point. The coexistence surface of the binary alloy model in the pa, Pb, T solid is 
shown schematically in Fig. 2. One is then led to ask whether the entire critical line 
shares a common behavior, or whether its character changes at some point. Although 
the former is clearly favored on the basis of universality, a careful examination of 
this question nevertheless appears worthwhile. 

As shown in Fig. 2, the Ising coexistence curve is the intersection of the critical 
surface with the plane p = Pa+Pb = 1- For p < 1 we can cross the critical line in the 
temperature or the density direction, maintaining all the while h = Pa~Pb = (pi is 
the chemical potential of species i) . Thus, while in the WRL there is no temperature 
per se, p = pa + Pb is a temperature-like variable, so that along the symmetry line 
h = 0, and in the vicinity of the critical point p c , the susceptibility should scale as 
X — (p — /U c )~ 7 , the order parameter, p A — Pb — (p — Pc) 13 , (for p > p c ), and so 
on PHI . (Note that if the chemical potentials are taken as independent variables, 
no "Fisher renormalization" of critical exponents is expected, as it would be if we 
worked with fixed densities fl33| .) 

Ref. f32| reports Monte Carlo simulations of the WRL in the grand canonical 



ensemble, using lattices of up to 160 2 sites in two dimensions and 64 3 sites in three 
dimensions. The algorithm employs three kinds of moves: "flips" (change of the 
state, A, B, or vacant, at a single site), "exchanges" (between any pair of sites in 
the system), and flips (A # B) of entire clusters (note that this is always possible 
in the WRL since a cluster of occupied sites is always surrounded by a border of 
vacant sites). The results for the reduced fourth cumulant, order parameter, and 
susceptibility are all consistent with Ising-like behavior p3 . 



9 



Following the WRL study, a series analysis of a closely related continuum model, 
with Gaussian Mayer /-functions, (i.e., the model introduced in Ref. |3l|), was re- 
ported by Lai and Fisher p^j . This study again provides good evidence for consis- 
tency with the Ising universality class at the critical point marking phase separation. 
More recently, the generalization of the WRL to q states, with infinite repulsion 
between unlike nearest-neighbor pairs has been found to exhibit critical behavior 
consistent with the g-state Potts model, as would be expected from symmetry con- 



siderations 35 



The appearance of Ising-like criticality in athermal models in fact goes back to 
studies of lattice gases with nearest-neighbor exclusion (NNE), done in the mid- 
1960's |36 H38ll . In this case there is only a single species of particle, and the only 



interaction is a hard-core repulsion assigning infinite energy to pairs of particles with 
separations < 1, in units of the lattice spacing. On a bipartite lattice in two or more 
dimensions, there is a critical density above which the particles begin to occupy one 
of the two sublattices preferentially, signalling a continuous transition to a state with 
antiferromagnetic order. (The NNE lattice gas is in fact closely related to the zero- 
temperature line of the Ising antiferromagnet.) Gaunt and Fisher analyzed series 
expansions for the NNE lattice gas on various two- and three-dimensional lattices, 
and found f3 ~ 1/8, suggesting, once again, Ising universality |36|,|37j]. It would be 
worthwhile applying modern series and simulation methods to the NNE models, to 
obtain more precise results on their critical behavior. 

It is perhaps worth observing that in all of the models mentioned in the present 
section (and, indeed, in the hard-sphere fluid), phase separation is driven exclusively 
by entropy maximization. Thus we have the apparently paradoxical conclusion that 
for large values of the chemical potential, the ordered phase has an entropy greater 
than or equal to that of the disordered phase. This shows that the naive identi- 
fication of entropy with 'disorder' is not always appropiate; interpreting increased 
entropy as greater freedom seems more apt. 



V. ISING WITHOUT EQUILIBRIUM: VOTERS, CATALYSIS, AND 

HYSTERESIS 

The present section concerns far-from-equilibrium models, defined by a Marko- 
vian dynamics rather than a Hamiltonian. Since the transition rates do not satisfy 
detailed balance with respect to any reasonable energy function (i.e., bounded be- 
low, and with a finite number of many-body terms), these models have no thermody- 
namic interpretation. Such systems are nevertheless widely studied using the tools 
of statistical physics, as models of, for example, populations, traffic, catalysis, and 
'self-organized" criticality |4(| . Many of these systems exhibit transitions between a 



fluctuation-free absorbing state (admitting no escape) and an active phase f4l|--f43 
Models possessing an absorbing state, and, associated with this, a non-negative 
order-parameter density, do not exhibit Ising symmetry. Instead, they fall generi- 
cally in the universality class of directed percolation, which plays a role analogous to 
the Ising class for absorbing-state phase transitions [p9| , ^4[ . 



The models to be discussed in this section do, however, observe Ising symmetry, 
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and belong to the Ising universality class. The fact that the Ising class cuts across 
the equilibrium/nonequilibrium boundary serves to illustrate that phase transitions 
in far-from-equilibrium systems are just as worthy of the name as their equilibrium 
counterparts, and do not, as is sometimes asserted, represent merely an "analogy" to 
"real" (i.e., equilibrium) phase transitions. Of course, one's outlook will depend on 
whether one adopts a thermodynamic or a mathematical definition of a phase tran- 
sition. My point is that experience with percolation, nonequilibrium models, and 
deterministic as well as stochastic cellular automata make it natural to regard, quite 
generally, any singular dependence of the properties of a system of many interacting 
units upon its control parameters as a phase transition. From this vantage, critical 
phenomena represent a particular class of singularities (attended, for example, by 
a diverging correlation length and relaxation time), with equilibrium critical points 
representing a particular, and not especially privileged, subset thereof. 



A. The Majority- Vote Model 

The prototypical example of a nonequilibrium stochastic process with local inter- 
actions, isotropy, and up-down symmetry is the majority- vote model MM. At each 



site of a lattice there is a "spin" variable <jj = ±1. A Markov process is defined as 
follows. At each step, a site i is chosen at random, and Oi — > which is taken to be 
equal to the majority of its nearest neighbors (i.e., sgn[J2jnni with probability 
p, and to the minority with probability q — 1 — p. (If there is no majority oi flips 
with probability 1/2.) Although the transition rates do not satisfy detailed balance, 
it is not hard to see that the parameter q plays a role analogous temperature in the 
equilibrium kinetic Ising model. (The dynamics of the majority vote model is akin 
to Glauber dynamics; there are no conserved quantities.) For q = 0, "voters" have 
no independene of opinion whatever, and slavishly follow the local majority. This 
q = limit defines the voter model, which has two absorbing states, all +1 or all 
-1, just as in the Ising model at T = 0. (The voter model belongs to a universality 
class different than Ising, the so-called compact directed percolation class p^j46|j .) 
For q > but small, the stationary state exhibits spontaneous "magnetization": 
m = (crj) = ±m , with m > 0. Above a critical q value (q c ~ 0.075 for the square 
lattice ^7j), the stationary state is disordered (m = 0). (No ordered state is found 
in one dimension, just as for the equilibrium Ising model.) 

Thus the phase diagram of the majority-vote model is qualitatively the same as 
for the standard Ising model. Confirmation that the critical behavior is also of the 



Ising kind comes from Monte Carlo simulations reported by de Oliveira ||47|| , which 
showed (for the two-dimensional case), that the exponents u, 7 and (3 are the same 
as those of the Ising model. Many nonequilibrium models have now been shown to 
exhibit Ising-like critical behavior, indicating that this universality class is equally 
robust, in or out of equilibrium. Examples include a family of generalized majority- 



vote models ||48|| , an anisotropic majority- vote model |49|], a two-state immunological 
model [ 50 1 , and various two-temperature Ising models PJ . In the latter example, a 
spin system evolves via two dynamical processes, for example, single-spin flips and 
nearest-neighbor exchanges, each having its own temperature (see Ch. 8 of Ref. 
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for a review). Ising-like short-time critical behavior has also been established for 



certain nonequilibrium models with up-down symmetry |52 . 

The principle underlying universality of critical behavior in these examples ap- 
pears, once again, to be symmetry While detailed renormalization group analyses 
are lacking, Grinstein et al. have argued that models obeying up-down symmetry 
and exhibiting a continuous phase transition, be they in equilibrium or not, will have 
a coarse-grained descrption of the same form as for a standard kinetic Ising model 
53] . We can see this as follows: if such a description exists, then in a Langevin-like 



equation for the order parameter density m(r, t), only odd powers of m may appear 
(up-down symmetry). Gradient terms are ruled out by isotropy, so the lowest-order 
derivative term will be V 2 m. Generically, the transition to an ordered state will be 
controlled by the coefficient a of the term oc m. (In mean-field theory this coefficient 
vanishes at the critical point, but in a full analysis the nonlinear terms will renor- 
malize a c to a nonzero value.) That is, the order-parameter density will obey (in 
the nonconserving case) a time-dependent Landau-Ginzburg equation of the form 

dm(r,t) _ 2 o . . ,_. 

— — — = V m — am — om + r/(r, t) , (7) 

where we have dropped higher-order terms (oc m 5 , etc.) since they are irrelevant to 
critical behavior. The noise term rj(r, t) is zero-mean, and Gaussian, with autocorre- 
lation {r)(r, t)r](r', t') = T5 d (r — r')5(t — t'). In equilibrium T is of course proportional 
to temperature. Out of equilibrium no such relation exists, but the theory again 
predicts an Ising-like transition; the deviation from equilibrium is found to be irrel- 
evant, near the upper critical dimension pBJ. The principle that systems with the 



same symmetries share a common critical behavior, whether in or out of equilibrium, 
has recently been extended to models with Potts-like symmetries ]5^J5B| . 

We have seen that a variety of nonequilibrium perturbations preserve the charac- 
ter of the transition in spin systems or lattice gases. It is worth noting that biasing 
the hopping rates to favor movement along a particular axis (e.g., in a lattice gas 
with attractive interactions and a Kawasaki-type exchange dynamics) breaks an 
essential symmetry of the Ising model by introducing a preferred direction. Such 
"driven diffusive systems" , as they have come to be called, exhibit critical behavior 
outside the Ising class. The precise nature of the transition remains controversial; 



for reviews see Refs. [43] and |56 . 



B. Catalysis Models 

Another large class of nonequilibrium lattice models arises in the study of het- 
erogeneous catalysis, typically on a metallic surface. One of the first such models to 
be studied in detail was introduced by Ziff, Gulari and Barshad (ZGB), to describe 
the reaction CO + 1/2 O2 — > CO on a platinum surface pT/| . These models typi- 
cally exhibit transitons to an absorbing state, which (when continuous) fall in the 
directed percolation class |5S[; they are reviewed in chapter 5 of Ref. P5[ . 

The ZGB model exhibits two phase transitions, one continuous, the other dis- 
continuous. Under a suitable perturbation, however, the latter can become a critical 
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point, which appears to belong to the Ising class. In the ZGB model, the catalytic 
surface is represented by a two dimensional lattice. The transition rates involve a 
single parameter, Y: the probability that a molecule arriving at the surface is CO. 
This molecule needs only a single vacant site to adsorb, as indicated by experimental 
studies; O2 requires a pair of vacant sites. After each adsorption event, the surface 
is immediately cleared of any CO-0 nearest-neighbor (NN) pairs, making the reac- 
tion rate for CO2 formation, in effect, infinite. The adsorption and reaction events 
comprising the dynamics go as follows. First choose the adsorbing species — CO 
with probability Y, 2 with probability 1 — Y — and a lattice site x (or, in the 
case of 2 , a NN pair, (x,y)), at random. If x is occupied (for 2 , if x and/or y 
are occupied), the adsorption attempt fails. If the newly-adsorbed molecule is CO, 
determine the set O(x) of NNs of x harboring an O atom. If this set is empty, the 
newly-arrived CO remains at x; otherwise it reacts, vacating x and one of the sites 
in 0(x), (chosen at random if 0(x) contains more than one site). If the newly- 
adsorbed molecule is 2 , construct C(x), the set of neighbors of x harboring CO, 
and similarly C(y). The O atom at x remains (reacts) if C(x) is empty (nonempty), 
and similarly for the atom at y. 

At a threshold value of Y, y 2 (— 0.5256 on the square lattice), the ZGB model is 
said to exhibit "CO-poisoning." That is, the stationary state changes abruptly from 
an active one (Y < y 2 ), with continuous production of C0 2 , to an absorbing state 
with all sites occupied by CO. (With all sites blocked, there is no way for 2 to 
adsorb.) The transition is sharply discontinuous. Experiments on surface-catalyzed 
CO + 1/2 O2 -> CO reactions do show a discontinuous transition between states of 
high and low reactivity as the partial pressure pco (analogous to Y in the model), is 
increased With increasing temperature, the transition softens, and at a certain 
temperature becomes continuous, about which the reaction rate becomes a smooth 
function of pco- This appears to be the result of thermally-activated, nonreactive 
desorption of CO. 

A similar sequence of alterations is observed if we include nonreactive desorption 
of CO at a certain rate, k, in the ZGB model. With CO desorption there is no 
longer a CO-poisoned state, but for small k the discontinuous transition between 
low and high CO coverages persists. Above a critical desorption rate k c ~ 0.0406, the 



coverages (and the C0 2 production rate) vary smoothly with Y |60H6~2| . Simulation 



results indicate that at the critical point of the CO transition, the critical exponent 
v = 1, and the reduced fourth cumulant u ~ 0.61, as would be expected for the 
two-dimensional Ising model W2 . 



C. Hysteresis in the Ising Model 

For our final example we come full circle to the Ising model, forced out of equilib- 
rium, this time, not through a bias or conflicting dynamics, but by a time- dependent 
external field h(t). Of principal interest is a periodically varying field, for example, 
h(t) sinusoidal or a square wave, with 

(/*) = ! J\(t)dt = 0, (8) 



1.3 



where T is the period of oscillation. Well below the critical temperature, and at suf- 
ficiently low frequencies, we expect to observe a hysteresis loop in the magnetization- 
field plane. As we raise the frequency, the system will at some point be unable to 
follow the rapid variations of the field, and the magnetization m(t) will remain close 
to one of its stationary values, i.e., +mo or — m . That is, for high frequencies we 
expect 

1 r T 

Q = — / m{t)dt ^ (9) 



r Jo 

whereas at low frequencies Q should be zero. 

For a sinusoidal field, the time-dependent mean-field equation describing this 
system predicts a continuous transition between a dynamically disordered (Q = 0) 
and a dynamically ordered (Q ^ 0) state at a critical frequency, which depends on 
the temperature and h = max[h(t)} |63| . The transition was confirmed in various 
Monte Carlo simulations, and may be relevant to experiments on thin magnetic films 
11,11 ■ 



Rikvold and coworkers suggested that the transition between dynamically or- 
dered and disordered states occurs when the period T becomes comparable to the 
metastable lifetime r(ho, T) of the state with the "wrong" sign of the magnetization 
(in a static field of magnitude ho), and presented extensive simulation results to 
support this proposal |56fl . These authors made a detailed numerical study of the 



critical behavior of Q in two dimensions; their results for exponent ratios and the 
fourth cumulant are consistent with Ising values. While the very general symmetry 
arguments outlined above again favor such a conclusion, finding a path from the 
microscopic dynamics of this nonstationary process to an equation of the form of 
Eq. (0) seems a particularly challenging theoretical task. 



VI. SUMMARY 

I have reviewed a number of examples in which critical behavior in the Ising 
universality class appears, despite the presence of long-range interactions, lack of a 
temperature scale, or lack of equilibrium. The Ornstein-Zernike approach that comes 
closest to reproducing Ising-like behavior was also reviewed. A general conclusion 
is that symmetries - of the order parameter and of the dynamics - are essential 
to determining critical behavior; the question of equilibrium is of minimal, if any, 
significance. On the other hand, we have seen that the effective, coarse-grained 
dynamics can look quite different from the microscopic interactions. In many cases 
there remains a considerable gap between these two levels of description. 
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Figure Captions 

Fig. 1. Best estimates for the location of the critical line and coexistence curve in 
the lattice RPM (simple cubic lattice). 

Fig. 2. Schematic coexistence surface of the binary alloy model. The Ising model 
coexistence curve is the intersection of the coexistence surface with the plane pa + 
p B — 1. The WRL coexistence curve is the intersection with the plane T = 0. 
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